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ABSTRACT 

Reliability  is  an  important  engineering  requirement  for  consistently  delivering  acceptable 
product  performance  through  time.  It  also  affects  the  scheduling  for  preventive  maintenance. 
Reliability  usually  degrades  with  time  increasing  therefore,  the  lifecycle  cost  due  to  more 
frequent  failures  which  result  in  increased  warranty  costs,  costly  repairs  and  loss  of  market  share. 
In  a  lifecycle  cost  based  design,  we  must  account  for  product  quality  and  preventive  maintenance 
using  time-dependent  reliability.  Quality  is  a  measure  of  our  confidence  that  the  product 
conforms  to  specifications  as  it  leaves  the  factory.  For  a  repairable  system,  preventive 
maintenance  is  scheduled  to  avoid  failures,  unnecessary  production  loss  and  safety  violations. 
This  article  proposes  a  methodology  to  obtain  the  optimal  scheduling  for  preventive  maintenance 
using  time-dependent  reliability  principles.  An  optimization  algorithm  maximizes  the  time  for 
preventive  maintenance  by  improving  the  system  reliability,  so  that  the  lifecycle  cost  stays  below 
a  specified  target.  The  lifecycle  cost  includes  a  production,  an  inspection,  and  an  expected 
variable  cost.  All  costs  depend  on  quality  and/or  reliability.  Preventive  maintenance  is  perfonned 
at  the  time  when  the  improved  reliability  falls  below  an  acceptable  reliability  target.  The 
methodology  also  identifies  the  most  critical  component(s),  or  failure  modes,  which  if  repaired, 
will  improve  the  system  reliability  the  most.  We  demonstrate  the  proposed  approach  using  an 
automotive  roller  clutch  example. 


1.  INTRODUCTION 

Customers  and  product  manufacturers  demand  continued  functionality  of  complex 
equipment  and  processes.  Degradation  of  material  properties  in  time  and  stochastic  operating 
conditions  such  as  loading,  result  in  a  lowered  resale  value  due  to  inconsistent  perfonnance  and 
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even  premature  failure,  increased  warranty  cost  and  customer  annoyance.  In  product  design,  it  is 
desirable  to  consider  the  attributes  of  survivability,  quality,  reliability,  maintainability,  and 
safety,  and  reduce  scrap,  rework,  and  maintenance  time.  This  can  be  achieved  by  selecting 
dimensions  and  materials  so  that  the  system  performance  satisfies  all  design  specifications  over 
the  planned  life-time.  An  important  factor  in  a  customer’s  buying  decision  is  the  total  lifecycle 
cost  which  can  be  affected  considerably  by  a  high  maintenance  cost. 

Quality  and  reliability  are  two  important  attributes  that  must  be  addressed 
simultaneously.  Quality  measures  the  conformance  to  design  specifications  as  the  product  leaves 
the  factory  to  begin  its  lifecycle.  In  this  work,  quality  is  quantified  by  the  probability  that  the 
performance  measures  meet  the  specifications  at  the  initial  time  of  its  lifecycle;  i.e.  right  after  it 
is  produced.  Reliability  is  defined  as  the  probability  that  the  system  will  perform  its  intended 
function  for  a  specified  interval  of  time,  under  stated  operating  and  environmental  conditions. 
Reliability  is,  therefore,  related  to  product  functionality  which  is  detennined  by  the  so-called 
“hard”  and  “soft”  failures  [1,  2].  In  a  hard  failure  the  system  loses  functionality  due  to  a 
complete  breakdown,  while  in  a  soft  failure  the  system  is  functional  but  one  or  more 
performance  measures  are  out  of  conformance. 

Maintenance  is  also  an  important  attribute.  It  involves  any  activity  that  can  restore  or 
recover  a  system  should  it  breaks  down  or  gets  out  of  order.  This  type  of  maintenance,  known  as 
repair  or  unscheduled  maintenance,  can  be  very  costly  and  should  be  avoided.  In  contrast,  a 
scheduled  maintenance  occurs  before  any  failure  or  malfunction  in  order  to  keep  the  system  in 
working  order  preventing  a  potentially  serious  damage.  Preventive  maintenance  is  therefore,  a 
schedule  of  planned  maintenance  actions  performed  in  an  attempt  to  avoid  failures,  unnecessary 
production  loss  and  safety  violations. 

In  preventive  maintenance,  when  the  maintenance  frequency  increases,  the  deterioration 
failure  rate  is  reduced  [3].  As  a  system  cannot  be  maintained  at  all  times,  a  way  is  needed  to 
decide  when  inspection/maintenance  is  needed.  To  avoid  unnecessary  maintenance,  which 
increases  the  cost  and  may  also  increase  the  risk  of  further  wear,  the  effectiveness  of  a  preventive 
maintenance  schedule  should  depend  on  quality,  reliability  and  a  lifecycle  cost  analysis.  Barlow 
and  Hunter  [4]  did  a  pioneering  work  on  optimum  preventive  maintenance  policies.  Since  then, 
many  works  have  been  done  along  this  line.  Survey  papers  such  as  in  [3-8]  cover  the  subject  of 
optimal  maintenance. 

In  our  previous  work  [9],  a  design  methodology  was  presented  based  on  the  product 
lifecycle  cost  considering  uncertainty  and  variability  through  time-dependent  reliability  in  order 
to  simultaneously  meet  the  needs  of  the  manufacturer  (low  production,  inspection  and  warranty 
costs)  and  the  customer  (customer  satisfaction  and  better  resale  value).  Although  design  under 
uncertainty,  or  reliability-based  design  (RBD),  has  been  extensively  studied  [10-14],  RBD 
considering  the  lifecycle  cost  has  not  attracted  much  attention  due  to  the  complexity  and  high 
computational  effort  in  calculating  the  time-dependent  reliability.  Only  a  limited  number  of 
studies  exist  which  mainly  use  simplifying  assumptions  in  the  calculation  of  time-dependent 
reliability  [1,2,  15-17]  or  use  lifecycle  cost  in  design  to  account  for  reliability,  maintenance  and 
warranty  [18-20]. 

The  intent  of  this  paper  is  to  increase  the  time  before  preventive  maintenance  is  needed 
using  design  optimization,  considering  time-dependent  phenomena  such  as  degradation  of 
certain  components.  Consider  a  system  with  multiple  responses  which  must  confonn  to  certain 
specification  limits  set  by  the  designer.  For  demonstration  purposes,  degradation  in  time  changes 
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the  statistics  of  the  inputs  (e.g.  operating  loads,  material  properties,  etc.)  and  the  resulting 
statistics  of  the  responses.  Fig.  1  shows  a  probabilistic  view  of  this  process  over  time. 


Input  Variables 


o 

a 

<D 

<D 

t-H 

Ph 


System  Responses  System  Reliability 

System 
Reliability 

Quality 

_ Reliability  (  r=  0) 
,  Acceptable 
Reliability 


tM  Time  (  r) 


f 

Time  for 
Maintenance 


Figure  1.  Uncertainty  propagation  for  time-dependent  system  reliability  (modified  from 

Savage  et  al  [1]) 

At  the  left  of  Fig.  1,  the  probability  density  functions  of  the  input  random  variables  Xt 
change  over  time  due  to  degradation.  They  are  therefore,  denoted  by  A((r)  to  indicate  that  they 
are  time-dependent.  The  designer  may  control  only  the  initial  probability  density  functions.  The 
system  response  distributions  also  change  over  time  due  to  the  changes  in  the  inputs.  It  is 
assumed  that  each  response  Z^r)  is  a  function  of  the  inputsX1(r),X2(r),...,X„(z').  The  system 

reliability  versus  time  is  shown  at  the  right  of  Fig.  1.  It  is  built  over  time  by  calculating  the 
system  time-dependent  reliability  or  probability  of  conformance.  As  a  special  case,  quality  can 
be  viewed  as  the  initial  reliability  at  r  =  0 .  If  an  acceptable  reliability  level  is  prescribed,  the  time 
to  maintenance  tu  can  be  determined.  The  efficient  estimation  of  the  time  dependent  reliability  is 
thus  the  key  to  determining  the  time  to  maintenance.  Maximizing  the  time  to  maintenance  for  a 
given  level  of  system  reliability  is  the  goal  of  this  research. 

To  date,  few  researchers  have  proposed  methods  to  estimate  the  system  reliability  by 
discretizing  time  and  using  time  invariant  limit-state  methods.  Using  a  different  approach,  Son 
and  Savage  [21]  tracked  time-variant,  limit-state  surfaces  in  the  standard  normal  space. 
Incremental  failure  probabilities  were  evaluated  and  summed  at  discrete  time  intervals  to 
estimate  the  mean  time  to  soft  failure  (MTTSF)  and  its  variance.  It  is  easier,  however,  for 
managers  and  engineers  to  understand  monetary  measures  rather  than  the  traditional  reliability 
measures.  Thus,  monetary  measures  are  commonly  used  as  optimization  objectives  or 
constraints.  The  costs  of  degrading  systems  can  be  used  to  detennine  the  means  and  tolerances. 

With  respect  to  performance  reliability,  there  have  been  a  few  research  activities  in  which 
the  present  value  (or  worth)  of  expected  quality  losses  are  used.  The  latter  are  measures  of  non¬ 
conformance  of  the  responses  with  respect  to  desired  specifications.  Chou  et  al.  [22,  23] 
augmented  present  value  with  the  production  cost  CP  to  form  a  total  cost,  and  then  perfonned 
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tolerance  allocation  by  formulating  an  optimization  problem.  Son  et  al.  [24]  used  a  probabilistic 
expected  loss  of  quality  and  extended  its  application  to  degrading  mechanistic  systems  with 
multiple  target/smaller/larger-is-best  performance  metrics.  The  probabilistic  measure  was 
optimized  using  aspects  of  Reliability-Based  Design  Optimization  (RBDO)  methods. 

In  our  previous  research  [9],  we  used  a  time-variant  reliability  analysis  and  loss  of  quality 
costs  [2]  in  the  design  of  multi-response  systems.  The  optimal  design  was  calculated  by 
minimizing  the  lifecycle  cost  while  satisfying  quality  and  reliability  constraints  throughout  the 
lifecycle  of  the  system.  The  confonnance  of  multiple  responses  was  treated  in  a  series-system 
fashion.  In  this  research,  we  maximize  the  first  time  for  preventive  maintenance  subject  to  an 
acceptable  quality,  time-dependent  reliability  and  cost  constraints.  An  automotive  roller  clutch 
example  illustrates  our  approach  and  provides  details  on  the  usefulness  and  practicality  of  the 
proposed  methodology  for  optimal  preventive  maintenance  scheduling. 

The  paper  is  organized  as  follows.  After  a  brief  introduction  to  time-dependent  reliability 
analysis  in  Section  2,  the  definition  of  lifecycle  cost  and  its  calculation  are  provided  in  Section  3. 
Section  4  presents  the  optimization  problem  for  maximizing  the  time  to  maintenance  and  Section 
5  uses  an  automotive  roller  clutch  example  to  demonstrate  the  benefits  of  the  proposed 
methodology.  Finally,  a  summary  and  conclusions  are  provided  in  Section  6. 

2.  RELIABILITY  ANALYSIS 

In  reliability  analysis,  we  estimate  the  probability  of  failure  due  to  variability  and/or 
uncertainty  using  a  probabilistic  quantification  of  the  variation/uncertainty  and  a  definition  of 
failure.  In  time-independent  (or  time-invariant)  reliability  analysis,  the  uncertainty  is  quantified 
using  random  variables.  A  time-independent  limit  state  function  g(x)  is  used,  where  the  vector 
X  represents  the  input  random  variables  with  a  joint  probability  density  function  /x(x)  .  The 
probability  of  failure  is  given  by 

P/=P(g(x)<0)=f...f  /s(x)A.  (1) 

Because  the  multi-dimensional  integral  of  Eq.  (1)  is  difficult  to  evaluate,  various  analytical  and 
simulation-based  methods  have  been  developed  [10-14,  25-30]. 

2.1  TIME-DEPENDENT  RELIABILITY  ESTIMATION 

In  many  engineering  systems,  the  system  response  (output)  depends  on  time  and  is 
described  by  a  random  process  g(d,X,r).  For  each  realization  of  the  input  deterministic  variables 
d  and  random  variables  X,  the  output  is  a  time-dependent  sample  function.  Denoting  the  time 
interval  of  interest  by  [ tmia ,  /rnax  J ,  the  probability  of  failure  within  this  time  interval  is  defined  as 

P} (d> X;  tmin,tmax)=p(3t  e  [tmin,tmax],  g(d,X,t)<  0),  (2a) 


or 
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Fc  (d,  X,  t)  =  P}  (d,  X;  0,  t)  =  P(3  r  e  [0,  t\  g(d,  X,  r )  <  0) 


(2b) 


where  Fc  denotes  what  is  known  as  cumulative  probability  of  failure  [30],  indicating  that  failure 
occurs  if  the  limit  state  function  becomes  negative  at  least  at  one  instance  within  the  time 
interval.  In  general,  the  definition  of  the  cumulative  probability  of  failure  in  Eq.  (2)  is  different 
from  the  instantaneous  probability  of  failure 

P;(d,X,t)  =  P(g(d,X,t)<0),  (3) 

which  is  calculated  in  a  quasi-steady  sense  by  fixing  time  t  and  replacing  the  random  process  by 
a  random  variable.  As  a  result,  it  does  not  account  for  possible  failures  before  timet .  The 
conceptual  difference  between  the  cumulative  and  the  instantaneous  probabilities  of  failure  is 
very  important. 

In  our  previous  work  [9],  we  gave  a  brief  literature  review  of  time-dependent  reliability 
evaluation  methods,  and  presented  a  composite  limit  state  concept  to  ease  the  estimation  of  the 
time-dependent  reliability  in  Eq.  (2).  The  composite  limit  state  in  the  time  period  [0,  t\  is  defined 
as 


gc(d,X,t)=  {3r  €  [0,4  g(d,X,r)<  0} »  U  (g(d,X,r/)<  0 


/=o 


(4) 


where  t,  e  M  for  r0  =  0  <  zx  <...<  t1  <...<  tl  =  t.  If  L  ->  oo  ,  the  approximation  in  Eq.  (4) 
becomes  an  equality. 

The  event  of  first  passage  in  a  given  time  interval  is  equivalent  to  the  event  of  the 
maximum  response  in  the  same  time  interval  exceeding  an  allowable  threshold  value.  The 
definition  of  the  cumulative  probability  of  failure  Fc{t)  is  thus  equivalent  to  the  probability  of 
first  passage.  If  the  maximum  response  in  the  time  period  [0,  t]  is 

ymax(d,X,t)  =  max  v(d,X,  r)  ,  (5) 

0  <T<t 

the  time-independent  “composite”  limit  state  for  a  fixed  time  interval  [0,  t]  is  defined  as 
gc(d,X,t)  =  l  -  vmax(d,X,t)/  vl  =  0  where  y‘  is  the  allowable  threshold  value.  The  composite 
failure  is  then  indicated  by 

gc(d,X,t)  =  1  - vmax(d,X,t)/vl  <0.  (6) 

In  this  work,  we  use  the  composite  limit  state  concept  in  the  illustrative  example  of 
Section  5.  We  calculate  the  most  probable  points  (MPP)  and  linearize  each  of  the  original  limit 
states  at  the  corresponding  MPP,  and  then  use  Ditlevesen’s  narrow  bounds  [31]  to  estimate  the 
cumulative  probability  of  failure.  References  [9,  31]  provide  more  details. 
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3.  COST  MODEL 


The  lifecycle  cost  C,  includes  the  production  cost  Cp ,  the  inspection  cost  C,  and  an 
expected  variable  cost  C£  [2,  32],  These  costs  are  a  function  of  quality  and  reliability.  The 
lifecycle  cost  is  expressed  as 

CL  (d,  X,  tf  ,r)=  CP (d,  X)  +  C,  (d,  X,  t0 )  +  CFV  (d,  X,  tf  ,r),  (7) 

where  d  is  the  vector  of  detenninistic  design  variables,  X  is  the  vector  of  random  design 
variables,  t  f  is  the  product  life,  and  r  is  the  interest  rate  to  translate  all  future  costs  to  a  present 

value.  Note  that  all  components  of  CL  depend  on  d  and  X.  The  inspection  cost  represents  the 
cost  to  inspect  or  fix  initial  defects.  The  expected  variable  cost  Cf£  is  a  function  of  the  time- 

dependent  reliability  which  is  used  to  estimate  the  expected  present  value  of  repairing  and/or 
replacing  the  failed  products  during  the  product  lifetime. 

The  calculation  of  CP  and  C7  costs  is  straightforward.  The  expected  variable  cost  is  due 
to  loss  of  quality  and  is  given  by  [2,  32] 

cv(tf)=\cf{T)fT{T)e~rTdT  (8) 

0 

where  cf{v)  is  the  cost  of  failure  at  time  r  ,  e~rT  is  the  discount  factor,  and  fT{r )  is  the 

probability  density  function  (PDF),  conditional  on  the  event  that  the  system  is  in  a  conforming 
state  at  the  initial  time  r  =  0.  The  latter  is  calculated  as  the  time  derivative  of  the  cumulative 
probability  of  the  failure  Fc(d,X,r);  i.e. 

/rW=<*tei.  ,9) 

dr 

Clearly,  Cy(tf )  depends  on  the  time-dependent  probability  of  failure,  which  is  expressed  by  its 
PDF  fT(r).  The  calculation  of  fT (r)  is  a  significant  challenge  and  was  addressed  in  [9]. 

4.  PROBLEM  STATEMENT 

It  is  desirable  to  increase  the  time  for  inspection  and  maintenance  of  a  repairable  system 
while  satisfying  reliability  and  total  cost  requirements.  In  order  to  deliver  “best  value,”  an 
optimization  must  determine  the  optimal  design  by  maximizing  the  time  to  maintenance  while 
keeping  the  system  reliability  above  a  specified  target  and  the  lifecycle  cost  below  an  acceptable 
level. 

Page  6  of  16  UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release. 


In  this  work,  we  consider  the  first  time  to  maintenance.  The  following  design 
optimization  problem  is  solved 


max  tM 
d,nx.ox^M 


(d,  p.x ,  ox ,  tM ,  r)  <  CL 

(10a) 

N'(d,X,t0)<l-i?'(t0) 

(10b) 

Fc{a,X,tM)<\-R‘{tM) 

(10c) 

(10d) 

Px,  ^  Px  ^  PX£/  >  <*XL  ^  °x  ^  °xc 

(lOe) 

where  the  time  to  preventive  maintenance,  tM  ,  is  maximized  by  varying  the  detenninistic  design 
variables  d,  and  the  means  px  and  standard  deviation  ox  of  the  random  variables  X.  In  Eq. 
(10a),  Cz(d,px,ox,tM,r)  comprises  the  mean  production  cost  Cp(d,px),  the  mean  inspection 
cost  C, (d,px,t0),  and  the  expected  variable  cost  Cy(d,]ix,tM,r)  up  to  tM  . 

Eq.  (10b)  ensures  that  an  acceptable  initial  quality  is  met.  F'(d,X,t0)  =  P'f{t0)  = 
p(g{a,x,t0)<  0)  is  the  instantaneous  CDF  for  quality  at  time  t0  =  0  and  R'(t)  is  the  system 
target  reliability  for  t  e  [t0,  tM\.  In  Eq.  (10c),  the  cumulative  probability  of  failure  F‘(d,  X,/)  at 
any  time  t  e  [t0,  tM\  must  be  less  that  an  acceptable  probability  of  failure  (l  -  R'  (tM  ))  at  tM .  Eqs 
(lOd)  and,  (lOe)  provide  lower  ( L )  and  upper  ( U)  bound  constraints  for  d,  px ,  and  ox  , 
respectively. 


5.  A  ROLLER  CLUTCH  EXAMPLE 

An  automotive  roller  clutch  is  used  to  demonstrate  how  the  proposed  formulation  of 
Section  4  can  be  used  in  design.  The  model  is  discussed  in  detail  in  [2,  33,  34].  The  clutch 
assembly  consists  of  a  hub,  a  cage,  and  four  rollers  made  of  1020  steel  as  well  as  four  springs 
(see  Fig.  2).  The  springs  push  the  rollers  against  the  cage  in  order  to  remain  in  contact  with  the 
cage  and  hub.  If  the  hub  turns  counter-clockwise,  the  rollers  bind  and  the  resulting  torque 
“locks”  the  cage  with  the  hub.  Otherwise,  the  rollers  slip  and  there  is  no  torque  conversion.  The 
main  features  of  the  design  include  the  hub  and  roller  diameters  D  and  d,  respectively,  the  cage 
outer  and  inner  diameters  B  and  A,  the  roller  length  L,  the  number  of  rollers  N,  and  the  material 
modulus  of  elasticity  E  and  Poisson’s  ratio  v . 
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Figure  2.  Roller  clutch  schematic 


A  proper  operation  is  ensured  by  three  performance  measures:  the  contact  angle  a  ,  the 
torque  capacity  r  ,  and  the  hoop  stress  ah .  The  contact  angle  ensures  proper  binding  and 
slipping.  If  the  value  of  the  angle  is  greater  than  the  upper  specification  limit  or  lesser  than  the 
lower  specification  limit,  the  clutch  does  not  work  correctly  and  it  must  be  reworked  or  scrapped. 
The  torque  capacity  is  provided  by  the  friction  between  the  hub,  roller  and  cage.  The  design  must 
provide  sufficient  torque  without  exceeding  material  compressive  strengths.  The  hoop  stress 
results  in  from  the  circumferential  roller  friction  force  on  the  cage.  It  must  be  kept  below  a 
maximum  to  prevent  fatigue  failure  in  the  pipe-like  cage. 

The  angle  a  ,  torque  capacity  r  ,  and  hoop  stress  ah  are  functions  of  the  design 
variables  [34],  and  are  given  by 


and 


a  -  cos  l(S) , 


(11a) 


/  \2 

t  =  NL\  — 


D2d 


V  ci  J 


4  (D  +  d) 


S 


-S2 


N 

(  Y 

Y  Dd  2 

S 

A 

(b2+a2' 

2  n 

vcU 

K(D  +  d)  j 

kb2-a2) 

(lib) 


(He) 


where 
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D  +  d 
A-d 


(Hd) 


is  a  geometric  ratio,  and 


(He) 


is  a  constant,  and  cr,  is  the  average  pennissible  contact  stress. 

The  design  variables  in  this  work  are  the  hub  diameter  D,  roller  diameter  d,  and  the 
inner  cage  diameter  A.  The  remaining  design  parameters  in  the  angle,  torque  and  stress  relations 
are  assumed  constant  with  L  =  80  mm,  B  =120  mm,  and  N  =  4.  The  properties  of  the  selected 
1020  steel  material  for  the  hub,  roller,  and  cage  are  v  =  0.29,  E  =  207  GPa,  and  crc  =  3,790 

MPa.  The  design  variables  are  assumed  to  be  statistically  independent  and  normally  distributed. 
Their  means  and  standard  deviations  are  assumed  to  be  within  the  following  bounds 


55.0973  <juD  <55.4973,  22.66  <  fid  <  23.06,  101.49  <  juA  <  101.89 , 

0.04  <  aD  <  0.08  ,  0.03  <od  <0.1,  0.07  <  a,  <  0.133  .  (12) 


The  three  design  variables  degrade  over  time.  The  surface  wear  on  the  clutch  parts 
causes  the  inner  diameter  A  of  the  cage  to  increase  over  time,  and  the  hub  and  roller  diameters  D 
and  d  to  decrease  over  time.  We  assume  the  wear  degradation  rate  is  constant  over  time  with  a 
wear  rate  of  k  =  2.5x10’4mm/year.  The  time  dependency  of  the  three  design  variables  is  thus 
expressed  by 

D(t)  =  D(l-kt),  d(t)  =  d(l- kt),  and^4(t)  =  A{\  +  ki).  (13) 

For  proper  operation,  the  angle  must  be  0. 1 1  radians  with  a  tolerance  of  ±  0.06 
radians.  Also,  the  torque  must  be  at  least  3000  Nm,  and  the  hoop  stress  must  be  at  most  400 
MPa.  The  following  four  constraints  express  therefore,  the  perfonnance  measures 

g{{p,d,A)=  0.05  -a  <  0  (14a) 

g2{p,d,A)—  a  -  0.17  <  0  (14b) 

g3  ( D ,  d.  A)  =  3000  -  t  <  0 

(14c) 

g4  (D,  d,  A)  =  ah  -400x1 06  <  0  . 


(14d) 

The  production  cost  is  defined  as  [34] 
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indicating  that  a  higher  precision  (smaller  tolerance)  increases  the  production  cost. 

We  assume  that  the  lifetime  is  tf  =  20  years  and  the  interest  rate  is  r  —  3%  .  Similarly 

to  [2],  we  define  the  inspection  cost  as  C7(px,ax)  =  csF'(jix,ax,t0  =  0)  where  the  cost  for 

4 

scrap  is  cs  =  $20  andF'(px,ox,t0  =  0)=  >  0)) .  The  expected  variable  cost 

i 

*/ 

is  given  by  Eq.  (8)  as  Cy  (t  f)  =  cf  j  fj{t)e  "  dt ,  where  the  rework  cost  is  assumed  constant  and 

o 

equal  to  =  $20 ,  and  the  PDF  f‘.  (/)  is  the  derivative  of  the  cumulative  CDF  F‘r  (/)  with 
respect  to  time  according  to  Eq.  (9).  The  expected  variable  cost  Cy(tf)  is  calculated  using  a 
trapezoidal  integration  scheme. 

In  this  example,  we  did  not  use  a  minimum  acceptable  initial  quality  as  Eq.  (10b)  of 
the  general  problem  formulation  indicates.  Also,  we  did  not  have  any  deterministic  design 
variables  d.  The  general  design  optimization  problem  of  Eq.  (10)  was  then  reduced  to 


max  tM 

l'\  -ht 

s.t.  CL  (px ,  <rx ,  tM ,  r)  <  C, 

(16a) 

Fc(px,ox,tM)<l-7?f(tM) 

(16b) 

Px,  ^  P  x  ^  PX(/  ,  °Xi  ^  °x  ^  Ox 

(16c) 

where  the  cumulative  failure  probability  of  the  system,  given  by 

f  4  h 

Fc(\ix,ax,tM)=P\  (J (g, (D, d ,  A)  <  0)  ,  is  time-dependent  because  D,  d,  and  A  depend  on  time 


(see  Eq.  13).  The  normally  distributed  random  design  vector  X  =  {D,  d,A}  has  a  mean 
px  =  {fJ-D > Md > Fa }  anc*  a  standard  deviation  ox  =  {aD,ad,aA}.  The  random  variables  D,  d,  and  A 
are  assumed  independent  because  the  kt  tenn  in  Eq.  (13)  is  very  small.  The  lower  and  upper 
bounds  for  px  and  ox  arc  provided  by  the  side  constraints  of  relation  (12),  and  the  initial  values 
for  px  and  ax  are  {55.30,  22.86,  101.69}  and  {0.06,  0.065,  0.1015},  respectively.  The  target 
reliability  at  the  time  for  maintenance  is  set  to  R'  {tM )  =  0.9  .  In  order  to  study  the  effect  of  the 
target  lifecycle  cost  C[ ,  we  solve  the  problem  for  C}=  ($18,  $19,  ...,  $28). 

For  C'L  =  $18,  an  optimal  design  for  the  maximum  time  to  maintenance  Zyi  =  7.01  years 
is  obtained  atjuD=  55.4946  mm,  jUd  =  22.7562  mm,  juA  =  101.49  mm,  aD=  0.08  mm,  crd=  0.064 
mm,  and  cr^  =0.111  mm.  The  hub  diameter  D  reaches  its  upper  bound  and  the  cage  internal 
diameter^  reaches  its  lower  bound.  The  constraint  (16b)  becomes  active  indicating  that  the 
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target  reliability  requirement  of  0.9  is  met.  The  four  limit  states  corresponding  to  the  four 
constraints  of  Eqs.  (14a)  through  (14d)  are  shown  in  Fig.  3  where  the  optimal  design  is  shown  in 
the  standard  nonnal  space  for  tQ  =  0  and  tM  =  7.01  years  by  projecting  the  four  limit  states  and 

the  corresponding  MPPs  onto  the  \ux  =  uD  w3  =  uA  }  space.  The  progression  in  time  of  each 

limit  state  is  shown  from  0  years  (thick  line)  to  7.01  years  with  an  increment  of  1  year.  The 
corresponding  MPPs  are: 

MPPl  =[-0.713,-1.09,0.940]  at ,  =  7.01  years, 

MPP2  =[0.901,1.387,-1.200]  aU=0years, 

MPP3  =[0.746, 1-151,  -1.000]  aU  =  0  years,  and 
MPP4  =[0.800,1.930,2.170]  aU  =  7.0iyears 

with  reliability  indices ^  =  1.607,  /?2=2.043,  /?3  =  1.698  and  /?4  =  3.012. 


Figure  3.  MPPs  and  limit  states  at  the  optimal  design  for  a  target  reliability  of  0.9  and  a  target 

cost  of  $18 

As  time  progresses,  the  limit  state  gx  =  0.05  -  a  <  0  expressing  the  minimum  contact 
angle  requirement,  becomes  more  critical  as  it  progressively  approaches  active  status,  and 
becomes  active  at  7.01  years  (tangent  to  the  target  reliability  circle  at  MPPi).  On  the  other  hand, 
the  limit  state  g2  —  a  —  0. 1 7  <  0  expressing  the  maximum  contact  angle  requirement,  is  more 
critical  (closer  to  active  status)  at  t  =  0  years,  attaining  its  MPPo  at  0  years.  The  limit  state  for  the 
minimum  torque  requirement  g3  =  3000  -  r  <  0  is  more  critical  at  t  =  0  years  and  attains  its 
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MPP3  at  0  years.  Finally,  the  limit  state  for  the  maximum  hoop  stress  requirement 
g4  =<Jh  -400xl06  <  0  becomes  more  critical  and  reaches  its  MPP4  at  7.01  years. 

Figs  4  through  6  show  the  optimization  results  for  a  target  lifecycle  cost  C[  from  $18 
to  $28  with  an  increment  of  $  1 .  The  target  system  reliability  at  the  time  for  maintenance  is  equal 
to  0.9.  For  the  initial  design  (dashed  red  line  in  Figs  4  and  5),  the  target  reliability  of  0.9  is 
reached  at  3.37  years  (Fig.  4)  and  the  corresponding  lifecycle  cost  is  approximately  $19  (see  Fig. 
5).  For  the  initial  design  therefore,  the  time  to  maintenance  is  ?om  =  3.37  years.  In  contrast,  for 
the  optimal  design  corresponding  to  a  target  cost  of  C[  =  $18  to  $28,  the  optimal  time  to 
maintenance  is  tM  =  7.01  to  12.77  years  which  is  much  higher  than  3.37  years.  This  is  achieved 
by  pushing  the  reliability  curve  at  the  optimal  design  to  the  right  (blue  lines  in  Fig.  4).  Fig.  4  also 
shows  that  for  a  given  time  to  maintenance,  the  reliability  of  the  optimal  design  increases  as  the 
target  lifecycle  cost  C[  increases. 

The  flat  part  of  the  blue  lines  in  Fig.  5  indicates  that  the  time  to  maintenance  can  be 
initially  increased  without  any  additional  cost.  This  happens  because  the  available  cost  is  not 
spent  on  the  non-critical  requirements  and  is  instead  directed  to  the  critical  requirements.  Thus, 
an  optimal  allocation  of  resources  is  achieved. 

Although  the  maximum  time  to  maintenance  can  be  increased  by  increasing  the  target 
cost,  it  is  ultimately  limited  by  other  constraints  as  shown  in  Figs.  4,  5  and  6.  Thus,  it  cannot  be 
further  increased  even  if  we  provide  more  resources.  Figs  5  and  6  clearly  show  that  the  constraint 
of  the  target  cost  remains  active  only  until  the  target  cost  is  raised  to  the  level  of  approximately 
$26.  For  higher  C[  values,  the  side  constraints  become  active  restricting  the  time  to  maintenance 
from  increasing,  even  if  more  resources  are  provided. 

This  practical  example  clearly  demonstrates  that  the  time  to  maintenance  can  be 
optimized  even  if  we  have  restricted  resources. 
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Figure  4.  Reliability  versus  time  and  optimal  time  to  maintenance  for  clutch  example 


Figure  5.  Cost  versus  time  and  optimal  time  to  maintenance  for  clutch  example 
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30 


28 


<M  (years)  ^  =  12.77  years 


Figure  6.  Cost  versus  time  to  maintenance  for  clutch  example 


6.  SUMMARY  AND  CONCLUSIONS 

Reliability  is  an  important  engineering  requirement  for  consistently  delivering  acceptable 
product  performance  through  time.  As  time  progresses,  the  product  may  fail  due  to  time 
phenomena  such  as  time-dependent  operating  conditions,  component  degradation,  etc.  The 
degradation  of  reliability  with  time  may  increase  the  lifecycle  cost  due  to  potential  warranty 
costs,  repairs  and  loss  of  market  share.  In  design  for  lifecycle  cost,  we  must  account  for  product 
quality,  time-dependent  reliability,  and  preventive  maintenance.  Quality  is  a  measure  of  our 
confidence  that  the  product  conforms  to  specifications  as  it  leaves  the  factory.  While  quality  is 
time-independent,  and  reliability  is  time-dependent.  Reliability  depends  on  1)  the  probability  that 
the  system  will  perform  its  intended  function  successfully  for  a  specified  interval  of  time  (no 
hard  failure),  and  2)  on  the  probability  that  the  system  response  will  not  exceed  an  objectionable 
by  the  customer  or  operator,  threshold  for  a  certain  time  period  (no  soft  failure). 

For  a  repairable  system,  preventive  maintenance  is  scheduled  to  avoid  failures, 
unnecessary  production  loss  and  safety  violations.  We  presented  a  methodology  in  this  article,  to 
increase  the  time  to  preventive  maintenance  while  satisfying  a  minimum  system  reliability  and  a 
maximum  allowable  lifecycle  cost.  This  was  demonstrated  using  an  automotive  roller  clutch 
example. 
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